How much random a random network is : a random matrix analysis 
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We analyze complex networks under random matrix theory framework. Particularly, we show that 
A3 statistic, which gives information about the long range correlations among eigenvalues, provides 
a qualitative measure of randomness in networks. As networks deviate from the regular structure, 
A3 follows random matrix prediction of linear behavior, in semi-logarithmic scale with the slope of 
I/tt^, for the longer scale. 
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Triggered by recently available data on large real world 
networks (e.g. structure of the Internet or molecular net- 
works in the living cell) combined with increasing com- 
puter power, an avalanche of qualitative research on net- 
work structure and dynamics are currently stimulating 
diverse scientific fields such as food web, nervous sys- 
tems, cellular metabolism, scientific collaborations, In- 
ternet, human language etc. P, H, Q- Results obtained 
so far span areas from the prevention of computer viruses 
to the stability and diversity of systems such as Inter- 
net, regulatory circuits of genome and ecosystems. These 
real world networks have several universal features, like 
small diameter, large clustering coefficient, scale-free de- 
gree distribution, assortative or disassortative mixing of 
the nodes, module structures [4|, etc. Irrespective of real 
world networks having one or more above mentioned fea- 
tures, one thing is common in all of them, and that is 
the existence of some amount of randomness or disorder 
in the connections structure. According to many recent 
studies, randomness in connections is one of the most 
important and desirable ingredient for the proper func- 
tionality or the efficient performance of systems having 
underlying network structure [U, [E B S Hi ■ For instance, 
information processing in brain is considered to be highly 
influenced by random connections among different mod- 
ular structure other examples of recent studies sug- 
gest the importance of randomness in various fields, such 
as evolution of language capacities of different species 
, evolution of cooperation in game theory Q , etc. On 
one hand these studies emphasize on the randomness in 
systems, on the other hand importance of structure or 
regularity is known for functional performance [3| ■ 

Seminal work of Watts and Strogatz shows that a very 
small amount of randomness has drastic impact on the 
diameter of networkjleading to the so called small-world 
(SW) phenomenon Q. Starting with a regular one di- 
mensional lattice, random rewiring of connections leads 
to the SW behavior for very small rewiring probability 
p. This simple rewiring scheme preserves the regular lat- 
tice structure yielding high clustering coefficients, while 



increases in the number of random connections among 
distant nodes resulting in the decrease of diameter of the 
network. Though this model network captures two most 
important characteristics of the real world networks, the 
whole structure of the network still remains a very reg- 
ular one-d lattice type, whereas real world networks are 
far more complex with a large number of random con- 
nections or seemingly random connections. The question 
arises whether one can identify or characterize the level 
of randomness in the complex networks ? There could be 
various possible ways to look into the problem. Recently 
L. da F. Costa has attempted to characterize random- 
ness by looking for regular patterns in networks [l^ . The 
purpose of present work is to have a qualitative measure 
of randomness in networks using eigenvalue fluctuation 
statistics of underlying adjacency matrix. 

Eigenvalues distributions have been studied profusely 
in the literature [111, They have certain signa- 

tures of the underlying network structures, like for com- 
plete random network it follows Wigner semi-circular 
law, for scale-free it follows triangular shape and for 
the small- world networks spectral distribution has multi- 
peaks structure ■ In this paper we are interested in 
comparing the randomness in the network structure and 
except for the above mentioned typical network struc- 
tures density distributions do not contain much informa- 
tions about the network structures, particularly it is not 
able to address the issue of "/low much randomness^^ . 

Our earlier contributions [H, [ij, [l^ have shown 
that the spacing distributions of complex random net- 
works follow universal behavior of random matrix theory 
(RMT). In this paper we characterize the level of random- 
ness in networks by using spectral rigidity test of RMT. 
Range for which A3 statistics follows random matrix pre- 
diction increases with the increase in random connections 
in the network. 

RMT was proposed by Wigner to explain the statisti- 
cal properties of nuclear spectra [l^ . Later this theory 
was successfully applied in the study of different complex 
systems including disordered systems, quantum chaotic 
systems, spectra of large complex atoms, etc [l^- More 
recently, RMT has been applied successfully to analyze 
time-series data of stock-market JAjlftL atmosphere 20] , 
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human EEG [2lj | , and many more lUISl. Nearest neigh- 
bor spacing distribution (NNSD) of the eigenvalues fol- 



lows two universal properties depending upon the under- 
lying correlations among the eigenvalues. For correlated 
eigenvalues, the NNSD follows Wigner-Dyson formula of 
Gaussian orthogonal ensemble (GOE) statistics of RMT; 
whereas, the NNSD follows Poisson statistics of RMT for 
uncorrelated eigenvalues. 

We denote the eigenvalues of the adjacency matrix of 
network by Ai, i = 1, . . . , N, where N is the size of the 
network and A^+i > A^V i. Eigenvalues are unfolded us- 
ing the technique described in [3, [3, [13 ■ Using the 
unfolded spectra {A^}, we calculate the nearest- neighbor 
spacings as 

= A,+i - A,; 
The NNSD P{s) for the case of GOE statistics 

P{s) = -sexpi^- — y (1) 

The NNSD reflects only local correlations among the 
eigenvalues. The spectral rigidity, measured by the As- 
statistic of RMT, provides information about the long- 
range correlations among the eigenvalues and is more 
sensitive test for RMT properties of the matrix under 
investigation [3 HI] ■ In the following, we describe the 
procedure to calculate this quantity. 

The Aa-statistic measures the least-square deviation 
of the spectral staircase function representing the cumu- 
lative density A^(A) from the best straight line fitting for 
a finite interval L of the spectrum, i.e., 

1 /■^+^ _ _ 9 _ 

A3(L;a;) = -min / [N{X)-aX-b] dX (2) 

L a,b 

where a and b are obtained from a least-square fit. Av- 
erage over several choices of x gives the spectral rigidity 
A3(L). For GOE, A3(i) depends logarithmically on L, 
i.e., 

A3(L)~^lnL. (3) 

Following we describe the method to generate networks 
with different randomness. Starting with one dimension 
ring lattice of N nodes, in which every node is connected 
to its k/2 nearest neighbors, we randomly rewire each 
connection of the lattice with the probability p such that 
self and multiple connections are excluded. Thus p — 
corresponds to a regular network, and p — I gives a com- 
pletely random network. The typical small- world behav- 
ior is observed around = 0.002 ^li]. Our earlier papers 
[H [Hi show that at this value of p = Pc NNSD follows 
universal GOE prediction of RMT. Universal GOE be- 
havior suggests that there exists a minimal amount of 
randomness in the networks yielding to the short range 
correlations among the corresponding eigenvalues. With 
the increase in p, obviously randomness increases in the 
network, but NNSD is not able to provide this informa- 
tion of enhancement of randomness. 
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FIG. 1: (Color online) Figure shows changes in spectral den- 
sity with the increase in random connections in the network. 
Initially, for (a)-(d) p = 0.002, 0.005, 0.01 and 0.02 : spec- 
tral density show typical small-world network like distribu- 
tion without having any specific analytical form with multi- 
ple peaks, (e)-(g) p — 0.05, 0.1 and 0.2 : approaching to- 
wards well known Wigner-Dyson semicircular distribution of 
random network, but having very long tail for larger posi- 
tive eigenvalues, (h) p = 1.0 : completely random network, 
showing exact Wigner-Dyson semicircular distribution. All 
networks have = 2000 nodes and an average degree A; = 20 
per node. Figures are plotted for the average over 10 random 
realizations of networks. 

As discussed in the introduction, NNSD only gives the 
information about the correlations among the neighbor- 
ing eigenvalues. It does not tell anything about the long 
range correlations among the eigenvalues. We probe long 
range correlations among the eigenvalues using A3 statis- 
tic, which tells that how closely the network follows ideal 
behavior of GOE of random matrices. For different values 
oi p > Pc, we construct several networks of size N = 2000 
and average degree fc = 20 With increase in the value of 
p, number of random rewired connections increases. We 
study the spectral rigidity of networks generated for var- 
ious p and then present the results for ensemble average 
of networks for each p value. 

We start with the network generated for p = pc = 
0.002, below this value of p spectral properties of net- 
works can not be modeled by GOE statistics of RMT 
(see our earlier papers [11,111). At p = Pc, there exists 
minimal amount of randomness sufficient to create corre- 
lations among the eigenvalues yielding GOE statistics of 



3 




50 100 150 50 100 150 50 100 150 50 100 150 

L L L L 



FIG. 2; (Color online) Figure shows changes in spectral behavior for different values of p. Network parameters remain same as 
for ll]). (a) p = 0.002 : on the verge of the small- world transition, NNSD follows RMT predictions, A3(Z/) statistic following 
RMT prediction only upto L ~ 30. (b) p = 0.01 : spectral density Hlb)) and NNSD show similar behavior as for the previous 
p value while A3(I/) statistic following RMT prediction for larger scale, upto L ~ 40. (c) p = 0.02 : ^z[L) statistic following 
RMT prediction upto L ~ 75. (d)p = 0.05 : for such a small value of p, A3(L) statistic following RMT prediction for very 
large scale (L ~ 135). For A'^ = 2000, we can not have meaningful A3 statistic for the larger L (see text). Hence for p > 0.05, 
p = 0.1 0.2 and 1.0: ^s,[L) statistic following RMT prediction for the scale L ~ 135. 



NNSD. Figure[T]plots spectral distribution of several net- 
works generated by Watts-Strogatz algorithm with differ- 
ent p > Pc values. As p increases, local regular structure 
destroys and random connections among nodes increase. 

Figure m plots NNSD distribution and A3 statistic of 
eigenvalues spectra of these networks. A3(L) is calcu- 
lated following Eq. ^ respectively. It can be seen from 
the plots in Figs. [Ha)-(d), that after p = Pc NNSD re- 
mains same with universal GOE statistics and we 
do not infer anything more about the randomness in the 
network. Though NNSD does not address the question 
of how much randomness, it contains a very important 
information; that various other aspects of RMT can now 
be applied to study these networks. 

Figures (e)-(h) plot A3 statistic of several networks 
with different p values. As it is seen from the figures A3 
statistic has consistent different behavior for different p 
values. Figs [2I^a)-(d) are plotted in the increasing order 
of the p values, with (a) being plotted at the transition 
to the small-world behavior. It can be seen from Fig- 
ures [1] (a) - (c) , that the density distribution of eigenval- 
ues for p values (0.002 — 0.02) are very similar with typical 
multi-peaks type of structure shown by small world net- 
works. One can not infer from these figures that which 



network has high random connections and which one has 
less. NNSD plots are also same for these networks. But 
interesting information comes out when we look at the 
spectral rigidity plots. 

Figure ^e) shows that the spectral rigidity follows 
GOE predictions upto range L 30. As p is increased, 
the scale for which spectral rigidity follows GOE pre- 
diction, in general, increases. Figure (H^f) is plotted for 
p = 0.01, for this value of p spectral distribution and 
NNSD remains same for such a small difference of ran- 
domness from the previous p values, but spectral rigidity 
now shows a bit different length scale {L ~ 45) for which 
it follows GOE predictions. Note that all the figures are 
plotted for the ensemble average of 10 networks gener- 
ated for each p values, so the statistics is the generic 
behavior of the networks with that randomness. For p 
taking value 0.02, Figure [2] (g) shows that there exists a 
considerable change in the length scale for which spec- 
tral rigidity follows GOE prediction upto L ^ 75. There 
is a consistent increase in the value of L with the p, in- 
dicating that A3 provides good measure of randomness 
in the networks. Figure MJ^) is plotted for p = 0.05, 
at this p network follows GOE prediction till very long 
range L ~ 135. It means that eigenvalues which are 
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L = 135 distance apart, are also correlated. Note that 
as one probes for the larger range correlations, the sta- 
tistical error in the calculation of A3 increases [l^l, and 
for the size of network {N = 2000) we are considering 
here it becomes difficult to provide a meaningful test of 
RMT using this statistics. For p > 0.05, though spectral 
densities change (Figure mf)-(h)) being semicircular at 
p ~ 1, the spectral rigidity plots (FigslJ) (f)-(h) remains 
same as for p ^ 0.05. According to the RMT interpreta- 
tion, eigenvalues of the network for p ~ 0.05 are as much 
correlated as the network for p ~ 1, and for p = 0.05 net- 
work has as much randomness as in a complete random 
network. 

We study complex networks under random matrix the- 
ory framework and use A3 statistics to provide a measure 
of deviation from the regular structure. Starting with 
a regular lattice, i.e. a network with complete regular 
structure, some connections are randomly rewired with 
probability p. Larger values of p indicate more deviation 
from regular structure in the network. In this manner 
several networks are generated with different deviations 
from the regular structure. According to the RMT : A3 
statistic, which measures long range correlations among 
the eigenvalues, provides insight about the randomness 
in the corresponding matrix. For a complete random ma- 
trix, eigenvalues are correlated till very long range. We 
use this result of eigenvalues correlation statistics to dis- 
tinguish the randomness among the networks, and show 
that A3 can be used as a qualitative measure of ran- 
domness. With the increase in the value of p, the length 



scale for which A3 follows GOE prediction of RMT also 
increases. In the semi-logarithmic plots, the slope of A3 
(Eq. [3]) matches exactly the GOE predictions of '--^ 1 /tt^ . 

Interestingly, for the networks generated withp ^ 0.05, 
A3 statistic already starts following GOE prediction for 
the scale as long as for the networks with p ^ 1. Ac- 
cording to the RMT, this tells that the eigenvalues of 
the networks generated with p ~ 0.05 are as much cor- 
related as for the networks generated with p ^ 1, which 
are complete random networks. We interprete this result 
as following: for p ~ 0.05 there is some kind of spreading 
of randomness over the network yielding the correlations 
among the eigenvalues as large as for the networks with 
P = 

To conclude, continuing our RMT analysis of complex 
networks, this paper uses A3 statistic to study the ran- 
domness in the networks. It is shown that RMT is partic- 
ularly suitable to study this important aspects of complex 
networks, thus on one hand enhancing the applicability 
of RMT, on the other hand opening a completely new 
domain to the complex networks studies. We hope that 
the complexity of the systems psj having network struc- 
tures can be addressed under the RMT framework, and 
thus helping to understand the dynamical behavior (26j 
and robustness of such networks better [131 . 
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Berlin) for encouraging discussions on the importance of 
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(currently at MPIPKS, Dresden) from random matrices 
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